g-Breathers in finite lattices: nonlinearity and weak disorder 
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Nonlinearity and disorder are the recognized ingredients of the lattice vibrational dynamics, the 
factors that could be diminished, but never excluded. We generalize the concept of g-breathers - 
periodic orbits in nonlinear lattices, exponentially localized in the reciprocal linear mode space - 
to the case of weak disorder, taking the Fermi-Pasta-Ulan chain as an example. We show, that 
these nonlinear vibrational modes remain exponentially localized near the central mode and stable, 
provided the disorder is sufficiently small. The instability threshold depends sensitively on a par- 
ticular realization of disorder and can be modified by specifically designed impurities. Basing on it, 
an approach to controlling the energy flow between the modes is proposed. The relevance to other 
model lattices and experimental miniature arrays is discussed. 

PACS numbers: 63.20.Pw, 63.20.Ry, 05.45.-a 



Nonlinearity and disorder are ubiquitous and unavoid- 
able features of discrete extended systems, the key play- 
ers in a wealth of fundamental dynamical and statistical 
physical phenomena such as thermalization, thermal con- 
ductivity, wave propagation, electron and phonon scat- 
tering. Studying lattice vibrational modes is vital to 
gain full understanding of these problems. Nonlinear- 
ity induces interaction between normal modes and en- 
ergy sharing if strong enough (known as the Fermi-Pasta- 
Ulam (FPU) problem) and time-periodic exponential 
localizations in direct space (discrete breathers) Lin- 
ear systems with disorder support another generic class 
of exponentially localized modes (Anderson modes) 13] • 
At the same time a satisfactory full understanding of the 
concurrent effect of nonlinearity and disorder is missing. 
This gap is being progressively filled for strongly disor- 
dered and weakly nonlinear lattices by intensive research 
on continuation of Anderson modes into nonlinear regime 
wavepacket spreading Q, light propagation in pho- 
tonic lattices 0, and Bose-Einstein condensate (BEC) 
localization in random optical potentials 

Little, however, is known on how the systems with pro- 
nounced nonlinearity and weak disorder behave. Remark- 
ably, it is a demand in a number of experimental and 
applicational contexts, beside a challenge from theory. 
Micro- and nano-electro-mechanical systems are rapidly 
developing components in microinstruments design (ul- 
trafast sensors, radio frequency filters, added mass sen- 
sors) [^. Their array structures offer broadband excita- 
tions, elastic waves, and effects of dispersion to be uti- 
lized Importantly, they are often suggested to 
operate in the nonlinear regime, while maturing technol- 
ogy reduces fabrication errors, hence diminishing spatial 
disorder. On the atomic scale, three-dimensional gold 
nano-cluster structures are found to be active and se- 
lective catalysts for a variety of chemical reactions, the 
surface vibrational modes being possible initiators 

One of the fundamental types of nonlinear oscillatory 
modes is q-breathers (QBs) - exact time-periodic solu- 
tions continued from linear modes and preserving ex- 



ponential localization in the linear normal mode coor- 
dinates. Originally proposed to explain the FPU para- 
dox (the energy locking in low-frequency modes of a 
weakly nonlinear chain, recurrencies, and size-dependent 
stochasticity thresholds) they have been discovered 
in two and three dimensional acoustic lattices, scaled to 
infinite systems, found in discrete nonlinear Schrodinger 
(DNLS) arrays, and quantum QBs were observed in the 
Bose-Hubbard chain [l3'|. Recently, QBs have been sug- 
gested as major actors in a BEC pulsating instability and 
a four- wave mixing process in a nonlinear crystal [lJ|. 

In this paper we extend the concept of QBs to ran- 
dom nonlinear media, exemplifying in the FPU chain. 
The cornerstones of our approach are continuation of 
QBs into non-zero 'frozen' disorder, taking a nonlinear 
localized solution as a seed, and statistical analysis of 
continued solutions. We show that QBs demonstrate the 
crossover from the exponential localization near the cen- 
tral mode to plateaus at a distance. The average sta- 
bility threshold in nonlinearity remains the same in the 
first order of approximation. In contrast, the standard 
deviation increases linearly with disorder, manifesting a 
high sensitivity on realizations. We analyze the effect of 
the harmonic in space inhomogeneities and discuss the 
energy fiow control by impurities design. 

The FPU-/? chain of N equal masses, coupled by 
springs with disorder in linear coefficients and quartic 
nonlinearity in potential, is described by the Hamilto- 
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where Xn{t) is the displacement of the n-th particle 
from its original position, Pn{t) its momentum, xq = 
XAT+i = 0, K„ e [—1/2,1/2] are random, uniformly 
distributed, and uncorrelated with = a\&n,m, 

a\ = 1/12 in our case. A canonical transformation 
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^n{t) = J'^iY, ^^(O sin defines the recipro- 

cal wave number space with N normal mode coordinates 
Qq{t), being solutions to the Hnear disorder-free case. 
The normal mode space is spanned by q and represents a 
chain similar to the situation in real space. The equations 
of motion read 
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are the normal mode frequen- 
l3| induce the se- 
lective nonlinear interaction between distant modes and 
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the all-to-all Hnear interaction due to disorder. 

Our methodology consists of two steps. Firstly, we 
take a known QB solution for non-zero nonlinearity 
[l^ . A particular realization of {k„} is chosen and 
d — D/^JN + 1 regarded as the disorder parameter. To- 
gether with the nonhnearity parameter v = j3/{N + 1), 
it is assumed to be small: u,d <^ 1. Then, an asymp- 
totic expansion in powers of {v, d\ is developed. Linear 
stability analysis is based on the constructed solution. 
Secondly, we address the statistical properties of the QB 
solution and instability threshold calculating respective 
averages and variances. 

Continuation of QBs to /3, L» 7^ from /3 7^ 0, L» 
employs the same technique as to /3 7^ 0, Z? = from 
(5 = D = Q For v^d << 1 and small amplitude ex- 
citations the g-oscillators get effectively decoupled, their 
harmonic energy Eq = ^ (^Q"^ + ^^Qq^ being almost con- 
served in time. Single g-oscillator excitations Eq ^ for 
q = qo only are trivial time-periodic and g-localized so- 
lutions for (3 — D ^ 0. 

For the disorder-free chain such periodic orbits can be 
continued into the nonlinear case at fixed total energy 
[1^ because the non-resonance condition nuj^ ^ ^q=Aqo 
(n being an integer) holds for any finite size [l5| and the 
Lyapunov theorem applies. Same ideas are expected 
to work for d ^ 0, as the spectrum remains non-resonant 
with the probability 1 Such continuation succeeded 
for all parameters we took. 

We continue QBs from nonlinear disorder-free solu- 
tions increasing D and keeping a particular random real- 
ization fixed. The total energy of the chain is E = 1 
in all examples, and 100 realizations of disorder are taken 
and used for averaging. Dependence of the average QB 
energy distribution on the level of disorder is reported 
in FigdJ We observe the characteristic QB exponentially 
localized profile on the almost fiat disorder-induced back- 
ground. The height of the plateau grows with D, grad- 



FIG. 1: The average mode energy distribution in QBs with 
go = 5,/3 = 0.01, Ai' = 32 under increasing disorder. Dashed 
lines are theoretical estimates (O. 



ually absorbing localized modes. For (3 = 0.01, go — 5, 
and N = 32 the second large mode is overcome near 
D = 0.01, and interpolation predicts the central one re- 
mains well above the background even for £> > 2, when 
some linear elasticity coefficients may become negative 
with a non-zero probability. 

Let us recall, that in case D = the QB solution 
Qq^{t) with a low-frequency seed mode number go can 
be written as an asymptotic ex pan sion in powers of the 
small nonlinearity parameter v |l2| . The energies of the 
modes go, 3go,. . . ,(2n -I- l)go,. . N read 
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and the frequency oj^^ = ujq^ (1 -I- %/4:vEq^^). Now we de- 
velop a perturbation theory to ([2]) in terms of the small 
disorder parameter d: Qq(t) — Qf'\t) + dQ[^\t) + . . ., 



its frequency being Cj = uj'^^' + du)'-'-' + . . ., substituting 
Qf\t) = Q^^{t) and w^") = u^^. In the first order 
approximation |[2]) becomes the equation of a forced os- 
cillator: Qg^^ + CJ^Qq^^ = —U!qUJq„Kq^qgQ''^g\ It folloWS 

that all modes get excited by disorder, their amplitude 
the bigger the closer its frequency to ujq^ : 



-Kq^qgAq„, q ^ qo, 
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the frequency being uj — ujqg{l + 9/4:iyEqg+d/2Kq,qg). As 
(Kq,qg) — 0, the first order corrections in d vanish for the 

averages 

(^q^) = {^^^^) = 0. Naturally, the variances 
are non-zero as the amplitude and frequency corrections 
vary depending on realization of {«;«}. The mode energy 
(averaged over period of a QB) approximately separates 
into nonlinearity and disorder-induced parts Eg « E^'" + 
°, where 



i^q^qo) ~~ of, .2 , .2 ^2 (^) 



Two limit cases are of particular interest: (i) q ^ qo, 
then « d^Eq„al/2 that gives a g-independent 

plateau energy (dashed lines in Fig[T]), and (ii) g = 
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FIG. 2: (a) The maximal and minimal absolute values of 
the eigenvalues of QBs with go = 6, A*' = 64 for two real- 
izations of {«:„}. For one realization the instability thresh- 
old in nonlinearity is increasing with D, for another - de- 
creasing, (b) Empty markers, dotted line: dependence of 
the variance of the instability threshold cr^* on the disor- 
der strength. Solid lines are theoretical estimates Filled 
markers: (/3*(D) - (/3*(0)). 
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Dctk < 2Eq„/{N + 1). Expressions (i) and (jS)) pre- 



diet the crossover between the exponential localization 
and the plateau at Qc ~ ^In ^{i^+i) I In-^ + 90 • Expres- 
sion (ii) suggests the 'small' Da^, ^ ^271^/ [N + lY and 
'large' Dctk ^ + 1) disorder criteria, the former 

implying that all QBs are localized and the latter that 
even QBs with go = 1 are delocalized. It also recovers 
the boundary oc \/N +1 between the localized (QBs) 
and delocalized in the q-space (but localized in the direct 
space) solutions (Anderson modes), that agrees with the 
previous results [l7|. Note, that the parameters taken in 
this paper correspond to the 'small' disorder. 

The linear stability of the continued periodic orbits is 
determined by linearizing the phase space flow around 
them and computing the eigenvalues 9i,i = 1, 27V of the 
corresponding symplectic Floquet matrix A QB is 
stable if \0i\ = 1, Vi. The maximal and minimal abso- 
lute values of 6i of QBs with go = 6, A'' = 64 for several 
increasing values of D and two different {ku} are plot- 
ted vs. /? in FiglSKa). Remarkably, while the instability 
threshold varies monotonously with D, it may not only 
decrease, but increase as well, depending on a particu- 
lar Moreover, stabilizing realizations are common, 
neatly balancing destabilizing ones. The observed devi- 
ation of the average instability threshold (/3*) from the 
disorder-free value /3q was much smaller that the variance 
(Figl2l(b)). The latter grows almost linearly in D, up to 
cr^j. K. 0.25/3(5, as seen for go = 6, TV = 64 (Figiajb); note 




FIG. 3: QB stability in case of spatially harmonic modula- 
tions: A'' = 32, D = 0.0025 and (a) p = 10, the central mode g 
is changed, (b) the central mode q = 8, the modulation wave 
number p is varied. Dashed lines are theoretical estimates. 



that for larger D the linear fit may become violated, due 
to the lower bound (3* > 0). 

The monotonous dependence of the instability thresh- 
old on D suggests that it is caused by the same resonance 
with the modes go ± 1 as in the disorder-free case. Let us 
explore the impact of disorder on this bifurcation. Lin- 
earizing equations of motion ([2]) around a QB solution 
Qq = Qq{t) + £.q{i), oue gets 



iq + ^liq = -SiyiVqEq^ COS^ {Ljt) ^ Cq^qg ^q^ ^pUJp^p 

P 

— dujq ''^ujpKp^g^p + 0{v'^, vd, d^) 
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The strongest instability is due to primary parametric 
resonance in ([H and involves a pair of the resonant modes 
g,p = go ± 1. Omitting non-resonant and 0{i'^,i^d,d'^) 
terms it is reduced to 



^q + Cj|(l + dKq^q)^q = -^l^UqUJpEq^ COS^ {'^t)S,p 
4 + w|(l + dKp^p)^p = -ivLOpUJqEqg COS^ [<^t)S^q 



(7) 



Thus, the disorder does not create new resonant terms, 
its impact being confined to the QB and resonant modes 
frequency shifts. The analysis analogous to yields 
the instability threshold /?*, its mean and variance: 



(/?*) = /3*, ap, = 2a,DVWTl/Eq,, 
where the disorder-free value is (3^ — 

AK = Kq^q - 2Kq„,q„ + Kp 

numerical results (Figl2]Jb)). 



(8) 



AK = Kq^q — 2Kqg,q„ + Kpp. It agrees well with the 
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One may ask now, which particular realizations fa- 
vor or disfavor stability? Furthermore, if some re- 
ceipts are disclosed, can they be used in controlling 
the energy flow in the mode space? The disorder- 
determined part of ([8]) can be rewritten as Ai^T — 

N+l 

~-7W+i i K„ cos ■^'^io(p-°-^^ sin^ ifj^+l) ■ It is linear 

n—l 

with respect to k„, thus we can represent it as a sum 
of spatial Fourier components, their contributions being 
additive. Thus, consider k„ = 0.5 cos ( ^^^^""'^^ +'p), 
where (p is the phase shift. It is natural to expect the min- 
imum of AK (and the maximal gain in stability), when 
p = 2qo, and it indeed yields AK — + 1 cos f, 

and the maximum P* = (3^ {l + D{N + i f /t:^) for 
(fi — 0. Immediately, a high sensitivity on ip is seen: 
the zero shift /3* = /^q for ip = ±7r/2; the minimum 
f3* = f3^{l- D{N + if/n'^) for <p = -tt. The ef- 
fect of p = 2(7o on adjacent QBs = go ± 1 is twice 
as small and reverse: for example, ii ip — Q then 
/?* = /3* (1 - D{N + l)V(27r2)). Remarkably, while for 
p = 2qo extremal shifts correspond to = 0, tt and zero 
ones to p — ±7r/2, for p = 2qQ ± 1 the zero shift ap- 
pears for (p — 0,7r, and the extrema for ip ~ ±tt/2: 
(3* = (3* (1 T 9,D{N + 1) V(37r3)) _ 

These results are illustrated in FiglH and show a good 
correspondence to the numerically determined QB sta- 
bility. That is, depending on the phase the spatially 
harmonic modulation of springs elasticities with the wave 
number p = 2qo, may significantly augment, weaken, or 
leave the stability intact (FiglSKa)). Modulations with 
p = 2qo ± 2 change the stability reversely and with 
twice a smaller amplitude for the same p, and those with 
p = 2qo ± 1 - just a bit weaker than 2qQ, but with a 7r/2 
shift in (p (FiglS^b)). Notably, modulations with other 
wave numbers have only a minor effect. Therefore, the 
spatial Foirier components with p S [2go — 2, 2qQ + 2] of 
{Kn} are decisive for the go-QB stability. 

These findings suggest a possibility of controlling the 
energy flow between modes. Indeed, by imposing a 
proper periodic modulation of the Hnear elasticity one 
can destabilize certain QB excitations and (i) promote 
equipartition or (ii) stabilize others, where the energy will 
be radiated; new QBs may also be subject to the same 
procedure to arrange the further energy flow. Experi- 
mentally, elasticity modulations could be achieved, for 
example, by laser heating, either as harmonic or spot im- 
purities, like it was designed to control discrete breathers 
location in cantilever arrays Q. 

In conclusion, we have demonstrated, that the concept 
of QBs can be successfully appHed to analyzing nonhnear 
vibrational modes in weakly disordered lattices. They es- 
sentially retain exponential localization and stability in 
the mode space, if the disorder is sufficiently small. We 
show, that the stability trend depends sensitively on a 
particular realization of disorder, and deliberately cre- 



ated inhomogeneities offer a promising technique of con- 
trolling the energy fiow between nonlinear modes. We 
expect that these ideas and methods to be applicable to 
a variety of nonlinear weakly disordered lattices - and 
we have already applied them to the DNLS chain (to 
be reported elsewhere) - including the contexts of a dif- 
ferent source of disorder (masses, nonlinearities), higher 
dimensions, and quantum arrays. The results on the non- 
linear modes sustainability, stability, and controUing are 
strongly expected to be in demand from experiments and 
applications. 

We thank S. Flach for stimulating and extremely valu- 
able discussions. 
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